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Abstract 

The isoscalar structure functions xF^ and F2 are measured as func- 
tions of X averaged over all permissible for the range of 6 to 28 GeV 
of incident neutrino (anti-neutrino) energy at the IHEP-JINR Neu- 
trino Detector. The QCD analysis of 0:^3 structure function provides 
= (411 ± 200) MeV under the assumption of QCD vahdity in the 
region of low Q^. The corresponding value of the strong interaction 
constant as{Mz) ~ ■123~^q'q\'^ agrees with the recent result of the 
CCFR collaboration and with the combined LEP/SLC result. 
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1 Introduction 



The data on deep-inelastic neutrino and anti-neutrino scattering in a wide 
region of momentum transfer provide a reliable basis for precise verification 
of QCD predictions [p. In this paper we present the measurements of the 
xF^ and F2 structure functions (SF) and the QCD analysis of in the kine- 
matic region of relatively small momentum transfer 0.55 < < 4.0 GeV'^. 
The value of the strong interaction constant as (Mz) is also evaluated and 
compared with the results of other experiments. 

2 Data Samples 

The analysis is based on data collected with three independent exposures 
of the IHEP-JINR Neutrino Detector Q to the wide-band neutrino and 
anti-neutrino beams |^] of the Serpukhov U-70 accelerator. The exposure 
to the anti-neutrino beam (l7^-exposure) was performed at the proton beam 
energy Ep = 70 GeV, whereas the two f^-exposures were carried out at 
Ep = 70 GeV and at Ep = 67 GeV. The energy of the resulting (v^) was 
in the range of 6 < -^1,(77) < 28 GeV. 

The experimental set-up and the selection criteria of charged current 
(CC) neutrino and anti-neutrino interactions are discussed in We re- 
stricted the range of measurements to W'^ > 1.7 GeV^ in order to reject 
quasi-elastic and resonance events and select mainly deep-inelastic neutrino 
and anti-neutrino interactions. The number of protons on target (p.o.t.) for 
each exposure, the selected number of CC and 77^ CC events and the 
mean values of for the three data samples are given in Table |l[ 



Table 1: Summary of the exposures. 



Beam 








Ep, GeV 


70 


70 


67 


Np.o.t. X 1017 


2.86 


1.05 


2.11 


Final statistics 


741 


2 139 


3 848 


{Q'), GeV2 


1.2 


2.3 
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3 Data Analysis 



The SF were measured as functions of x averaged over all permissible for 
the energy range 6 < E^^-p-j < 28 GeV. The events were binned in intervals 
of X, and the values of xF^ and F2 were calculated in these intervals. 

The number of interactions, n'^ , and 17^ interactions, n'^, in a given 
bin of X is a linear combination of the average values {F2} and {xF^} of the 
respective SF in this bin (we assume invariance under the charge conjuga- 
tion): 

= a^-{F2} - 6^-{xF3}, 

The subscripts 1 and 2 correspond to the f^-exposures at Ep = 70 GeV and 
Ep = 67 GeV respectively. The quantities a'^^'^^ and b'^^'^^ are the integrals 
( "flux integrals" ) of products of the differential neutrino (anti-neutrino) flux 
i^E) and the known factors depending on the scaling variables x and 
y, as given by the standard form of the differential cross-section for deep- 
inelastic Vfj_ (l7^)-scattering off an isoscalar target: 

vr 

= N- y{l- ^)E(l)^^''\E)dxdydE. 

vr J 2 

Here N is the number of nucleons in the fiducial volume of the detector 
and the parameter R = {F2 — 2xFi ) j 2xFi measures the violation of the 
Callan- Gross relation [^. 

The number n'^^'^^ of neutrino (anti-neutrino) interactions in a given bin 
of x was obtained from the measured number of neutrino (anti-neutrino) 
events in this bin corrected for acceptance, smearing effects arising from 
Fermi motion and measurement uncertainties, radiative effects (following 
the prescription given by De Riijula et al. Q) and target non-isoscalarity 
(assuming d^/uy = 0.5) Q. To determine appropriate correction factors, 
the Monte Carlo simulation of the experimental set-up has been carried out 
using the CATAS program We used the Buras and Gaemers (BEBC) 
parameterization Q for quark distributions. The charm quark content of the 
nucleon was assumed to be zero. The kinematic suppression of d ^ c and 
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s ^ c transitions was taken into account assuming slow rescaling [10| and 
the charm and strange quark masses of nic = 1.25 GeV and nig = 0.25 GeV 
respectively. The Fermi motion of nucleons was simulated according to jll| J . 
The details of the Monte Carlo simulation are described in 12 1. 

The number of interactions in a given bin of x is subjected to kinematic 
constraints imposed by the cuts on the muon momentum (p^ > 1 GeV/c 0]), 
on the neutrino (anti-neutrino) energy (6 < E^^-p^ < 28 GeV) and on the 
invariant mass square of the hadronic system {W^ > 1.7 GeV^). These 
constraints were taken into account in the calculation of the flux integrals 
by appropriate modification of the volume of integration. 

The measured values of F2 and xF^ structure functions are given in 
Table and in Fig. 1. 



Table 2: The isoscalar structure functions F2 and xF^ obtained with the 
assumption of i? = 0. The difference AF2 between the values of F2 obtained 
with i? = 0.1 and with 72 = is also presented. The bin edges are at x = 
0.02, 0.1, 0.2, 0.3, 0.4, 0.5 and 0.65. The shown systematic errors do not 
include the normalization error of 4% for F2 and 11% for xF^ originating 



from the uncertainties in the and i/^ flux prediction [13|. 



(x) 




F2 


stat 


syst 


AF2 


XF3 


stat 


syst 


0.052 


0.55 


1.169 


.039 


0.047 


0.023 


0.445 


.458 


0.062 


0.148 


1.4 


1.097 


.036 


0.022 


0.022 


0.583 


.087 


0.017 


0.248 


2.2 


0.894 


.032 


0.018 


0.019 


0.622 


.075 


0.019 


0.346 


2.9 


0.576 


.028 


0.017 


0.013 


0.556 


.109 


0.011 


0.447 


3.4 


0.390 


.025 


0.012 


0.009 


0.336 


.070 


0.007 


0.563 


4.0 


0.182 


.017 


0.004 


0.004 


0.177 


.117 


0.005 



The systematic errors presented in Table ^ come from the uncertainties 
in the knowledge of neutrino flux and cross-sections @ , imperfect detector 
calibration and uncertainties in the correction factors due to the choice of 
the input quark distributions. The overall normalization error, originating 
from the uncertainties in the and flux prediction |]l^] , was estimated to 
be 4% for F2 and 11% for xF^. The correction factor uncertainties were eval- 
uated by repeating the calculation of the SF using the Field-Feynman [0] 
and GRV |15] parameterizations of quark distributions. 
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The obtained experimental data on the xF^ were then compared with the 
QCD prediction for Q^-evolution by the Jacobi polynomials method in the 
next-to-leading order (NLO) QCD approximation |jl6|, 17, |l^. Performing 
the QCD analysis of the xF^ SF, we do not discuss here the problem of 
validity of application of perturbative QCD predictions for the kinematic 
region of low and do not take into account nuclear effects, heavy quarks 
threshold effects and higher order QCD corrections. 

In order to take into account the target mass corrections, the Nachtmann 
moments |19] of F3 are expanded in powers of M'^^^JQ'^. Retaining only 



the terms of the order of M'^^^JQ'^, one could obtain: 



M3(iV,Q' 



m: 



QCD 



.f2 



+ 



Here M^^^{N,Q'^) are the Melhn moments of xF^, : 



QCD 



X 



N-2 



xF3{x, Q"^) dx, 



N = 2,3,... 



(2) 



The Q^-evolution of Mellin moments is defined by QCD [20, 21 1 and is 
presented here for the non-singlet case for simplicity: 



m: 



QCD 



l3o 



as {Ql 
as (Q2 



dN 



3 

(0)NS 



2/3o, 



11 - -n. 



M?^^{N,Ql), N = 2,3, 



(3) 



Here as{Q'^) is the strong interaction constant, "y^^^^^ are the non-singlet 
leading order anomalous dimensions and is the number of flavours. The 
factor Hn {Qo,Q^) contains all next-to-leading order QCD corrections 



21, 22 . 
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The unknown coefficients M^'^'^ (NjQq) in (^) could be parameterized 
as the Melhn moments of some function: 

M^^^{N,Ql) = x^-^Ax\l-xydx, (4) 

JO 

iV = 2,3, ... 

where the constants A, b and c should be determined from the fit to the data. 
Having defined the moments (Q) - (^) and following the method discussed 
in lit, 17], we can write the xF^ SF in the form: 



max 

n=0 

n 
j=0 

where Q^'^{x) are the Jacobi polynomials and dj{a,j3) are the coefficients 
of the expansion of 0"'^(x) in powers of x: 

n 

@^'^ix) = J2c]{a,P)xr 

j=0 

The accuracy of the SF approximation better than 1% is achieved for 



^max = 9 in a wide region of the parameters a and /3 ||T7 ] . 



The higher-twist (HT) contribution is also taken into account: 
xFs{x,Q')=xF^''''{x,Q') + ^ 



where h{x) = 0.166 - 3.746 • x + 9.922 • x^ - 6.730 • is chosen by an 
interpolation of the NLO result for the HT contribution from [^]. This 
shape of h{x) is in a good agreement with the theoretical prediction of 



and with the result of 26 1 obtained for a higher kinematic region. 

Using nine Mellin moments and taking into account target mass correc- 
tions, we have determined four parameters - A, b, c and the QCD parameter 
Ajg^ (Table ^). In order to decrease the number of free parameters we have 
fixed the value of parameter A using the Gross-Llewellyn Smith sum rule 

- dependence: Sgls = 3 (1 - 0]. The fit was performed using the 

MINUIT program |28|. Three sources of errors - statistical, systematic and 
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Table 3: The results of the NLO QCD fit to the xF^ SF data for n/ = 4, 
= 3 GeV2, Nmax = 9, a = 0.7, (3 = 3.0. 





0.22 


A 


10.4 (fixed) 


b 


0.86 ± 0.14 


c 


3.83 ± 0.61 


MS 


(411 ± 200) MeV 


as{Mz) 


n 1 9q + U.UiU 

u.iza_o.oi3 



normalization - were summed up in quadrature. The errors for the free pa- 
rameters corresponding to the 70% confidence level were obtained using the 
procedure described in . A relatively good accuracy of the measurements 
of was achieved due to a high sensitivity of the QCD evolution equa- 
tions to the variations of A^g^ in the low region (0.55 < < 4.0 GeV^ 
in our case). 

The value of as{Mz) corresponding to the measured value of A-p^ 
was calculated from the so-called "matching relation" ||3^ and found to 
be as{Mz) = 0.123l[5:0}0. 

4 Discussion of the results 

We have compared our results with the measurements performed by other 
experiments. The comparison led to the following comments: 



The parameter A^ = (411 ± 200) MeV of the fit to the xF^ SF 



MS 

data is in agreement with the NLO analyses with HT contribution of 

MS 



the CCFR xFs data: A^ = (381 ±53 {stat)± 17 {HT)) MeV Q and 



= (428 ± 158 (exp)) MeV |2 

The value of obtained from the NLO analysis of the xFs SF 

provides the value of the strong interaction constant at the point 
of Z boson mass of as{Mz) = 0.123^00^3 which is in agreement 
with the result of the analysis of the CCFR data aa{Mz) = 0.119 ± 



0.002 (exp) ± 0.004 (f/ieory) H and with the combined LEP/SLC 
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result as{Mz) = 0.124 it 0.0043 |33|. Our current measurement is 
higher by one standard deviation than the value as{Mz) = 0.113 ± 
0.003 (exp) ±0.004 (theory) |3^, obtained in the F2 structure function 
analysis of the BCDMS and SLAG data on fiN and eA^ deep-inelastic 
scattering. 



5 Conclusion 

We have presented the measurements of the structure functions F2 and xF^ 
in the kinematic range 0.02 <x< 0.65 and 0.55 < < 4.0 GeV^, obtained 
from the inclusive deep-inelastic u^^ and 17^ scattering data collected at the 
IHEP-JINR Neutrino Detector. The NLO QCD analysis with HT contribu- 
tions of xF^ under the assumption of QCD validity in the region of low 
provides A^i^ = (411 it 200) MeV; the corresponding value of the strong 
interaction constant is as{Mz) = 0.123]'^o;gJ3. 
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Figure 1: The measured x-dependence of the isoscalar structure functions F2{x) and 
xF-^{x). The statistical and systematic errors are added in quadrature (the normahzation 
errors of 4% for F2 and 11% for xF^ are not shown). 
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